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We theoretically study the nuclear spin induced decoherence of a quantum dot in Si that is 
confined at a SiGe interface. We calculate decoherence time dependence on ^^Ge in the barrier 
layer to evaluate the importance of Ge as well as Si enrichment for long decoherence times. We 
use atomistic tight-binding modeling for an accurate account of the electron wavefunction which 
is particularly important for determining the contact hyperfine interactions with the Ge nuclear 
spins. We find decoherence times due to Ge spins at natural concentrations to be milliseconds. 
This suggests SiGe/Si quantum dot devices employing enriched Si will require enriched Ge as well 
in order to benefit from long coherence times. We provide a comparison of T2 times for various 
fractions of nonzero spin isotopes of Si and Ge. 
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I. INTRODUCTION 

An electron spin in silicon is a promising candidate 
for a quantum bit. Sources of decoherence such as 
spectral diffusion of nuclear spins can be nearly elim- 
inated through isotopic enrichment and other decoher- 
ence mechanisms such as spin-orbit coupling are small. 
Quantum dots offer a way to control and couple spins 
electrically, demonstrated in GaAs, and have shown fast 
single qubit gate operations. Significant progress towards 
achieving coherent control of spins in single and dou- 
ble quantum dots in Si/SiGe have been suggested re- 
cently [1-5 . Long decoherence times in silicon have been 
achieved [6], in large part, through enrichment of the Si 
with the zero nuclear spin isotope ^^Si. In SiGe/Si quan- 
tum dots, the wavefunction is in close proximity with a 
distribution of germanium nuclei in the SiGe barrier re- 
gion. Approximately 7.8 percent of the isotopes in natu- 
rally occurring pure germanium carry spin (i.e., ^^Ge). In 
this article we use a tight-binding method, NEM03D [7 , 
to calculate the spatial extent of the Si/SiGe quantum 
dot wavefunction and spin-bath cluster methods to calcu- 
late the decoherence behavior resulting from the presence 
of these non-zero, 9/2, nuclear spins. 



II. ELECTRON WAVEFUNCTION 
CALCULATION 

The spatial extent of an electron localized in a quan- 
tum dot is calculated using an atomistic tight binding 
(TB) approach embedded in the Nanoelectronic Model- 
ing Tool (NEM03D) [7 , which solves the single elec- 
tron eigen spectra for multi-million atom systems. We 
adopt this detailed atomistic approach to model a realis- 
tic quantum dot as the TB method in general captures a 
variety of atomic scale effects. In particular, this quan- 
titative description is needed to accurately obtain the 
wavefunction amplitudes on the Ge atoms, responsible 
for the contact hyperfine interaction [8] . A good descrip- 
tion of the spatial distribution of the wavefunction is also 



necessary to compute the dipolar interactions precisely 

We simulate a quantum dot confined vertically in a 
SiGe/Si quantum well and laterally by a harmonic po- 
tential, as described by the electronic (el) Hamiltonian, 

H^^ = H^^ + a{x - xof + a{y - y^f + eF,z (1) 

where Hq is the TB Hamiltonian of the Si-SiGe lattice. 
The 2nd and 3rd terms represent a 2D parabolic quan- 
tum dot of curvature a = 3.137 x 10~^ eV/nm^, which 
translates to a characteristic ground state energy E = 0.5 
meV and a dot radius of R = 2b nm (i.e. the distance 
at which the wavefunction amplitude falls off as 1/e of 
the peak value). The last term in Eq. ([T]) is due to a 
variable vertical E-field F^, which simulates the effect of 
a top enhancement gate [9]. A harmonic potential rep- 
resents a reasonably good estimate of the confinement 
for electrostatically gated structures [O]. T While mi- 
nor numerical variations of the results presented here are 
expected for much smaller dot sizes, in which the local 
alloy disorder has greater impact from sample to sample, 
the same trends and order of magnitude estimates will 
remain valid. 

The Hamiltonian of the host materials Hq is formu- 
lated semi-empirically with the 10 band sp^d^s* nearest 
neighbor model [7] to represent the full bandstructure of 
the hetero-structure. The TB model parameters are well- 
established in literature [12], and validated in a variety 
of earlier works [SJ [11] . The model of the SiGe region is 
purely atomistic [12], and no virtual crystal approxima- 
tions are involved. 

The vertical structure is a 10 nm thick, tensile strained 
silicon layer sandwiched between fully relaxed SiGe bar- 
rier layers of thickness 41 nm each, above and below the 
strained Si layer respectively, as shown in Fig. 1(a). The 
equilibrium lattice constants and strain in the two layers 
are also calculated through minimization of the strain 
energy calculated using the valence-force field Keating 
model ^]. The size of the lattice is chosen so that 
the SiGe, to a good approximation, is fully relaxed [Tij . 
For strain relaxation, simulations of 22 million atoms 
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were performed, whereas electronic structure simulations 
could be done on a reduced domain of 3.5 million atoms 
incorporating only a part of the SiGe regions near the 
strained Si region. The band offsets between Sii-^cGca^ 
and Si obtained from atomistic tight-binding simulations 
compare well with experimentally known values [15]. 

The full TB Hamiltonian, including the electrostatic 
potential given in Eq. ([T]) , is solved by a parallel Lanczos 
algorithm to obtain a number of the lowest energy states 
of the quantum dot. The ground state wavefunction of 
the QD is of relevance in this work. Since Si has a six- 
fold degenerate conduction band, the orbital states are 
multiplied by six in general. However, presence of strain, 
vertical E-field, and the SiGe/Si interface break the con- 
duction band valley degeneracy of the system [16l [17], 
introducing an energy gap between the ground state and 
the excited states. The splitting between the lowest two 
valleys increases with the vertical E-field [I6l [18] . Our 
consideration in this Article is for a spin qubit encod- 
ing and we assume that the excited valley states are kept 
much further than ksT above the two energy levels of the 
qubit states. We calculate valley splitting in our model 
to range from 0.1 to 1.2 meV in correspondence with our 
E-field range of 1 — 20 MV/m. This splitting is ample 
for double quantum dot qubits encoded in the spin zero 
subspace of singlet and triplet states [19 , which is the 
predominant system of consideration for silicon quantum 
dot based quantum computing. 
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FIG. 1: a) The simulated device structure showing the SiGe 
and Si regions for both strain and electronic structure calcu- 
lations, b) The expectation value of the wavefunction in the 
SiGe region as a function of E-field and Ge concentration x in 
Sii-ccGccc. c) Wavefunction probability density plot for the yz 
plane showing confinement at the interface, d) Wavefunction 
probability density in the xy plane near the interface inside 
SiGe showing variations due to alloy disorder, c) and d) are 
plotted for Fz = 20 MV/m and a Ge concentration x = 0.3. 

The penetration and overlap of the electron wavefunc- 
tion in the Sii-a^Ge^^ barrier is calculated as a function of 
both the Ge concentration x in the alloy and the vertical 



field, as shown in Fig.[T](b). The overlap is defined as the 
expectation value of the wavefunction in the Sii-a^Ge^; re- 
gion. The overlap magnitude changes 5 fold for the entire 
range of common experimental parameters examined. As 
expected, the more the Ge concentration, the higher is 
the interfacial barrier, and less the wavefunction overlap. 
Applying a higher E-field pulls the wavefunction more 
into the SiGe barrier, and increases the overlap. In Fig. 
2(c) and 2(d), we show two different cross-sections of the 
wavefunction, in the xz plane, and in the xy plane near 
the interface, respectively. The random fluctuations in 
Fig. 2(d) is due to the alloy disorder in SiGe. The dipole 
and hyperfine coupling to the nuclear spin bath in the 
Ge are calculated using the electron wavefunction enve- 
lope from the TB solution for each corresponding E-field 
strength and Ge concentration. 



III. SPIN BATH CLUSTER EXPANSION 
CALCULATIONS 

We examine the central spin decoherence of the quan- 
tum dot single electron spin in the silicon quantum well 
among remotely located ^^Ge nuclear spins in the SiGe 
barrier. The quantum dot spin experiences fluctuations 
of its energy splitting due to spectral diffusion (SD). To 
calculate the evolution of the central spin, we use a clus- 
ter expansion technique [20 , which has been shown previ- 
ously to agree extremely well with decoherence measure- 
ments of bulk silicon systems in the presence of varying 
concentrations of background ^^Si [2T| [22]. 

Here we independently consider the effects of various 
concentrations of either ^^Ge or ^^Si nuclear spins. We 
consider isotropic (contact) hyperfine and anisotropic hy- 
perfine (AHF) interactions of the quantum dot electron 
spin with the nuclei. We also consider dipolar interac- 
tions between the nuclear spins which are responsible for 
the dynamics of the nuclear spin bath that causes SD. 
We neglect quadrupole interactions which could be sig- 
nificant for ^^Ge; these interactions are unknown to our 
knowledge and its effects decrease with increasing mag- 
netic field. Taking the limit of a large magnetic field such 
that Zeeman energies are conserved and only fiip-fiop dy- 
namics of bath spins are allowed, the effective spin (sp) 
Hamiltonian is 
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written in atomic units; a factor of fi/(47reo) is implied 
for the Hamiltonian. We use S for denoted electron spin 
operators and In for denoted nuclear spin operators. For 
the gyromagnetic ratios, we use 7e = {g = 2)iib/^ = 
1.76 X 10^^(T s)-\ 7n = 0.936 x 10^(T s)-^ for ^^Ge, 
and 7n = -5.319 x 10^(T s)-^ for ^^Si. The hyperfine 
interaction, An-, consists of an isotropic Fermi contact hy- 
perfine term proportional to the probability of the elec- 
tron being on top of the nth nucleus and an AHF term 
that is the dipolar interaction averaged over the rest of 
the electron wavefunction. This AHF term is negligible 
compared with the contact hyperfine for nuclei within the 
quantum dot region but dominates far from the quan- 
tum dot. For |^(Rn)P of the contact hyperfine interac- 
tion, we must multiply the envelope of the wavefunction 
from the TB solution by a factor r] = ||ii(Rn)|p, where 
li(Rn) is the periodic Bloch function. For the ^^Si, we 
use T] = 186 [23 . For the ^^Ge, we use an r] that is ten 
times larger as a rough estimate [24 . The probability 
density is higher at ^^Ge sites than at ^^Si since Ge has 
more charge. The ^^Si are spin 1/2 nuclei and the ^^Ge 
are spin 9/2 nuclei. 

In the large magnetic field limit, the AHF contribu- 
tions from the Bloch part of the wavefunctions have a 
negligible effect. Relaxing this limit, these AHF contri- 
butions cause B-field dependent electron spin echo enve- 
lope modulations (ESEEM) [25] arising from S^I^ terms 
that can be significant in Si [26] . These Bloch-dependent 
AHF interactions strengths are calculated to be no more 
than 3% of the corresponding (per nucleus) contact hy- 
perfine f23^ strengths. At a B-field of 100 mT, we es- 
timate the maximum ESEEM amplitude to be at most 
3% in our models with natural Ge (0.1% due to natural 
Si with its smaller 77). This amplitude maximum scales 
with nuclear concentration and inversely with the B-field 
squared. Any other B-field dependence from keeping the 
hyperfine and dipolar interaction terms that do not con- 
serve Zeeman energy are negligible above 10 mT accord- 
ing to our cluster calculations. In our subsequent results, 
we neglect B-field dependent effects. 

The Hahn spin echo is used as a standard decoher- 
ence measure, T2 commonly reported as the time that the 
Hahn spin echo reaches exp(— 1) times the peak value. 
We compute Hahn spin echo decays and associated T2 
using a cluster expansion technique [20 . The Hahn spin 
echo consists of initializing the quantum dot electron spin 
to an equal superposition of up and down, freely evolv- 
ing for a time r, applying a spin-flipping tt rotation, then 
evolving again for time r. We treat the refocusing pulse, 
the midway tt rotation, as instantaneous and ideal. The 
echo decay is parameterized by the total evolution time 
t = 2r. The cluster expansion works particularly well 
where the interactions among bath spins are weak rel- 
ative to the interactions with the central spin, as is the 
case since 7^ <C 7e. We are therefore justified in using the 
lowest order of the cluster expansion, the pair approxima- 
tion, which accounts for the effects of flip-flopping among 
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FIG. 2: Hahn echo T2 times for various Ge concentrations 
and electric fields calculated for a SiGe/Si quantum dot elec- 
tron due to a naturally occurring abundance of ^^Ge and no 
^^Si. The applied magnetic field is taken to be perpendicular 
to the lateral direction. 



all pairs. This approximation is expressed as 



U = e-^^'er^^e-^^ff", (7) 

(aif/t|;)(t|f7|t))^^ (9) 

U((^\^krn\i){t\Un,rn\t)] 



Un,m 

VE{t = 2r) 



where VE{t) is the spin echo decay. We use heuristics to 
find the pairs with significant contributions and ignore 
the rest in order to make the calculation more efficient. 

Our theoretical Hahn echo results for various SiGe/Si 
quantum dot scenarios with decoherence induced by nat- 
ural Ge are shown in Fig. [2j We study the effects of dif- 
ferent Ge concentrations and different magnitudes of an 
applied electric field. The electric field pulls the quantum 
dot deeper into the SiGe layer so that T2 times decrease. 
An increase in concentration of the decoherence-inducing 
Ge causes T2 to decrease but this is mitigated to some 
extent by the corresponding reduction in the overlap of 
the quantum dot wavefunction into the SiGe [Fig. [l] (b)]. 
The Hahn echo decay due to natural Ge is approximately 
of the form exp(— t^). This behavior was theoretically 
predicted [20 , 27J in GaAs quantum dots and experimen- 
tally confirmed [28 . It is understood as the lowest order 
contribution in an expansion in time of the log of the 
echo. The decay is produced by a large number of small 
contributions from many different bath spin pairs, each 
well- approximated by this lowest order in time. This 
tends to occur in quantum dots that envelope a large 
number of bath spins but not in donors where there is 
a frozen core effect (nuclear spins near the center of a 
donor cannot flip-flop because of the difference in their 
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ppm '-"Ge or ^^Si 

FIG. 3: T2 and T2 for a SiGe/Si quantum dot electron as 
a function of the fraction of the nonzero spin isotope of Ge 
(^^Ge) in the barrier layer assuming spinless Si. Overlaid is 
the T2 and T2 of the quantum dot spin as a function of the 
fraction of the nonzero spin isotope of Si (^^Si) assuming spin- 
less Ge. The ^^Si results are indicated with thick black curves 
and the rest are ^^Ge results. T2 curves are indicated with 
star symbols. The magnetic field direction is perpendicular 
to the lateral direction for the solid curves and parallel to the 
lateral direction for the dotted curves. Colors (color online) 
indicate different electric field strengths. 

Knight shifts) [20]. Here, the bath spins are on the fringe 
of the quantum dot, but the effect is the same. 

^^Si enrichment has been used to attain long coherence 
times [6II2T]. It is also possible, in principle, to simply 
remove the decoherence-inducing Ge spins via isotopic 
enrichment. In Fig. [3j we examine T2 and times as a 
function of various fractions of ^^Ge or ^^Si nuclei. T2* is 
the exp (—1) decay time of a freely evolving ensemble of 
similar quantum dots with different random spatial and 
state configurations of the spin bath (the inhomogeneous 
broadening effect). For coherent control, gate operations 
should be fast relative to T2*. T2 is relevant for quantum 
memory. 

For a given ^^Si concentration it would be sensible to 
enrich the ^^Ge to attain a comparable T2. We observe 
that the effects of the electric field strength and magnetic 
field direction upon T2 are significant at higher concen- 
trations of ^^Ge, but are not so relevant at low concen- 
trations. In a sparse bath, such effects are minor. The 
magnetic field direction has little effect upon T2*. At 
moderate to high concentrations isotropic contact hyper- 
fine interactions dominate T2 and at low concentrations 
the system is almost rotationally invariant. The electric 
field, however, does impact T2* significantly since it has 
a strong impact upon the contact hyperfine interactions. 
At an electric field of 20 MV/m, the ^^Ge nearly ap- 



proaches that of ^^Si as a function of their respective iso- 
topic fractions. Although the overlap into SiGe at this 
field and x = 0.2 Ge concentration is only about 10%, 
this is counter-balanced by the order of magnitude in- 
crease of ?7 for the charge concentration at ^^Ge nuclear 
sites. 

The low concentration regime has a different charac- 
ter than the natural isotopic concentration regime. In 
the low concentration limit, the geometric details of the 
problem become irrelevent and the effective Hamiltonian 
scales as (all dipolar interactions) so that time 

scales inversely with spin concentration. This is denoted 
by the dashed line of Fig. [3j The generic problem in this 
regime is well-studied in Refs. [22l [29l From these stud- 
ies we note that there is substantial sample to sample 
variation such that the T2 of a given qubit may differ 
substantially from the ensemble average T2 that we re- 
port. 



IV. SUMMARY 

Decoherence of quantum dot electron spin in Si con- 
fined by SiGe barriers containing '^^Ge are calculated. 
The spatial extent of the wave function is calculated us- 
ing a tight binding approach, NEM03D, and the time 
dependence of the spin decoherence is calculated using a 
cluster expansion approach. The details of the SiGe pen- 
etration of the wavefunction are important for determin- 
ing the isotropic contact hyperfind and anisotropic hy- 
perfine interactions with the ^^Ge. Naturally abundant 
concentrations of ^^Ge are predicted to limit decoher- 
ence times, T2, to the order of milliseconds, many orders 
of magnitude below that experimentally demonstrated in 
enriched silicon (^^Si). The decoherence time is predicted 
to increase rapidly through enrichment with zero nuclear 
spin Ge isotopes and reach decoherence times longer than 
that demonstrated in silicon (about 1 second) at ^^Ge 
concentrations below 200 ppm. That is, other mecha- 
nisms in Si will likely limit the decoherence time at these 
enrichment levels. 
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